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Rational approximations to \J~2 and other 
algebraic numbers revisited 

par Paul M VOUTIER 

Resume. Dans cet article, nous ameliorons des mesures effec- 
tives d'irrationalite pour certains nombres de la forme \fn en 
utilisant des approximations obtenues a partir de fonctions hy- 
pergeometriques. Ces resultats sont tres proche du mieux que 
peut donner cette methode. Nous obtenons ces resultats grace a 
des informations arithmetiques tres precises sur les denominateurs 
des coefficients de ces fonctions hypergeometriques. 

Des ameliorations de bornes pour 9(k,l;x) et ip(k,l;x) (k = 
1, 3, 4, 6) sont aussi presentes. 

Abstract. In this paper, we establish improved effective irra- 
tionality measures for certain numbers of the form ^/n, using 
approximations obtained from hypergeometric functions. These 
results are very close to the best possible using this method. We 
are able to obtain these results by determining very precise arith- 
metic information about the denominators of the coefficients of 
these hypergeometric functions. 

Improved bounds for 6(k,l;x) and ip(k,l;x) for k — 1,3,4,6 
are also presented. 



1. Introduction 



In this article, we shall consider some refinements of a method due to 
Alan Baker [TJ [2] for obtaining effective irrationality measures for certain 
algebraic numbers of the form z m l n . As an example, he showed that for 
any integers p and q, with q ^ 0, 

10~ 6 



>l/3 _ V 



> 



12.955 ' 



There are two 

{-^m,n,r(-2 ; )} 



r=0 



infi- 
and 



This method has its basis in the work of Thue. 
nite families of hypergeometric polynomials in Q[z] 

{Y mi n,r(z)}^° =0 , such that {Xm,n,r[ Z ) I X m,n,r{ Z )}T=Q is a sec l uence of g°°d 

approximations to z m l n . Under certain conditions on z, these approxima- 
tions are good enough to enable us to establish an effective irrationality 
measure for z m l n which is better than the Liouville measure. 
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Since it is easy to obtain sharp estimates for the other quantities involved, 
the most important consideration in applying this method is the size of the 
denominators of these hyper geometric polynomials. 

Chudnovsky [5] improved on Baker's results by showing that, if p is a suf- 
ficiently large prime divisor of the least common denominator of X mnr {z) 
and Y m ^ n ^ r {z), then p must lie in certain congruence classes mod n and 
certain subintervals of [1, nr]. 

In the case of z m l n = 2 1//3 , he was able to show that for any e > there 
exists a positive integer qo(e) such that 

1 



2 l/3 



> 



12.4297.. 



\Q\ 

for all integers p and q with \q\ > qo(e). Moreover, since his estimates 
for the relevant quantities are asymptotically correct, this exponent is the 
best that one can obtain from this hyper geometric method although "off- 
diagonal" or the method of "ameliorating factors" (a la Hata) still might 
yield improvements. 

Shortly after this work, Easton [6] obtained explicit versions for the cube 
roots of various positive integers. For 2 1 ' 3 , he showed that 

2.2 • 1(T 8 



2 i/3 



P 



> 



12.795 



for all integers p and q with q ^ 0. 

It is the purpose of this paper to establish effective irrationality measures 
which come quite close to Chudnovsky's. In the particular case of 2 1 / 3 , 

0.25 



2 i/3 



V 



> 



12.4325 



for all integers p and q with q ^ 0. 

This paper was initially written and circulated in 1996. Independently, 
Bennett [3] obtained a result, which in the cubic case, is slightly weaker 
than the theorem stated here. E.g., for 2 1 / 3 , he showed that 

0.25 



2 i/3 _ P 



> 



12.45 



for all integers p and q with q ^ 0. 

In fact, this subject has been the topic of even more work. As part of his 
Ph.D. Thesis (see [8]), Heimonen has also obtained effective irrationality 
measures for numbers of the form y/a/b, as well as of the form log(a/fe). 
His results are not as sharp as those of the author, but they are still sub- 
stantially better than Easton's. 

The general method used in each of these three papers is essentially 
the same. However, there are substantial differences in the presentations 
due to the fact that the approach of Bennett and Heimonen shows more 
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apparently the role that Pade approximations play in this area, while the 
author deals explicitly with hypergeometric polynomials. 

Actually, the referee has pointed out that other work in this area has 
been done, producing results not much weaker than our own. And this 
work preceeded the results of Bennett, Heimonen and the author. We 
are referring to the work of Nikishin [12] and, especially, Korobov [9]. In 
particular, in 1990, Korobov showed that 



for all natural numbers p and q with g / 1,4. The reader looking for a 
more accessible reference to these works is referred to pp. 38-39]. 

The main differences between this version of the paper and the previous 
version are Theorem 12 . 31 and improvements in computer hardware. This has 
resulted in replacing 0.93 with 0.911 in the exponents on e in the expressions 
for E and Q in Theorem 12.11 (which requires a larger value of c\ ) , along 
with the consequent improvements to Corollary 12.21 including new results 
for tyiT and ^57. 




The main incentive for publication of this paper now is completeness. 
Several articles have since appeared in the literature (e.g., [TT] and [15] ) 
which depend on results in this article. Furthermore, the lemmas in this 
article, which are either new or sharpen results currently in the literature, 
are important in forthcoming articles by the author and others. They are 
accompanied by an analysis showing that they are best-possible or else 
what the best-possible results should be. And lastly, the main theorem 
itself, along with its corollary, is an improvement on the present results in 
the literature. 

We structure this paper as follows. Section 2 contains the statements of 
our results. In Section 3, we state and prove the arithmetic results that 
we obtain for the coefficients of the hypergeometric polynomials. Section 4 
is devoted to the proof of Theorem 12.3^ as this theorem will be required 
in Section 5, where we obtain the analytic bounds that we will require for 
the proof of Theorem 12.11 Section 6 contains the diophantine lemma that 
allows us to obtain an effective irrationality measure from a sequence of 
good approximations. At this point, we have all the pieces that we need 
to prove Theorem 12.11 which is done in Section 7. Finally, Corollary 12.21 is 
proven in Section 8. 

Finally, I'd like to thank Gary Walsh for his encouragement and moti- 
vation to resume my work in this area. Also, Clemens Heuberger deserves 
my thanks for his careful reading of an earlier version of this paper and 
accompanying suggestions. And, of course, I thank the referee for their 
time and effort as well as their suggestions for improvements. 



\/2 - p/q > q 



-2.5 
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2. Results 

Theorem 2.1. Let a and b be integers satisfying < b < a. Define c±,d,E 
and k by 

*/3 /(a -b), 
if 3 || (a — b) and 
otherwise, 

2 




E 



e -0.911 3 d( a l/2_ 6 l/2 



log j e 0.911 3 -d( a l/2 + fc l/2)2j 



log£ 



and 



ci 



10 40(«+1) 



a. 



IfE> \ then 



(1) 



(a/b^-p/q 



> 



cm 



K+l 



for all integers p and q with q ^ 0. 



Remark. c\ grows quite rapidly as the absolute values of the arguments of 
the exponential functions in the definition of E approach their best possible 
value of vr^/3/6 = 0.9068 .... 

In the earlier version of this paper with 0.911 replaced by 0.93, we could 
have taken c\ = 10 7 ^ K+1 ^a. It is feasible to prove Theorem 12.11 with 0.911 
replaced by 0.91, but then we would have to take c\ = 10 86 *- K+1 ^a. 

The rate of growth is even more rapid as we continue to approach 0.9068. 
For example, with 0.907, c x > 10 2400 ( K+1 )a. 

As an application of Theorem l2.lt we give effective irrationality measures 
for all numbers of the form y/n where n is a cube-free rational integer 
with 2 < n < 100 and for which the hyper geometric method yields an 
improvement over the Liouville bound. 

Corollary 2.2. For the values of n given in Table 1, we have 

C-2 



\/n — p/q\ > 



\K+1 ' 



for all integers p and q with q ^ where C2 and k are the values corre- 
sponding to n in Table 1. 



Remark. If a be an irrational element of Q (\/n), then we can write 



a 



«3 \fn + «4 ' 
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n C2 K 


n C2 K 


n C2 K 


2 0.25 1.4325 


25 0.07 1.7567 


60 0.08 1.5670 


3 0.37 1.6974 


26 0.03 1.4860 


61 0.06 1.5193 


4 0.41 1.4325 


28 0.03 1.4813 


62 0.04 1.4646 


5 0.29 1.7567 


30 0.10 1.6689 


63 0.02 1.3943 


6 0.01 1.3216 


31 0.14 1.9288 


65 0.02 1.3929 


7 0.08 1.6717 


36 0.08 1.3216 


66 0.04 1.4610 


9 0.08 1.6974 


37 0.01 1.2472 


67 0.06 1.5125 


10 0.15 1.4157 


39 0.08 1.1848 


68 0.08 1.5562 


11 0.22 1.8725 


41 0.41 1.9956 


70 0.12 1.6314 


12 0.28 1.9099 


42 0.12 1.4186 


76 0.08 1.5154 


13 0.35 1.8266 


43 0.01 1.2890 


78 0.03 1.5729 


15 0.19 1.4964 


44 0.21 1.8164 


83 0.09 1.6898 


17 0.01 1.1996 


49 0.13 1.6717 


84 0.37 1.8797 


18 0.37 1.9099 


50 0.11 1.1962 


90 0.09 1.3751 


19 0.02 1.2718 


52 0.26 1.8901 


91 0.009 1.2583 


20 0.009 1.1961 


57 0.15 1.9825 


98 0.38 1.4813 


22 0.07 1.2764 


58 0.12 1.6526 


100 0.35 1.4158 



Table 1. Results for 



where 01,02,03,04 S Z with 0104 — 0203 7^ 0. In this way, we can use 
Corollary 12.21 to obtain effective irrationality measures for any such a (see 
Section 8 of [5]). 

These values of a and b were found from the convergents p/q in the 
continued-fraction expansion of ^fn by setting a/b to be either (p/q) 3 /n or 
its reciprocal, whichever is greater than one. For each cube-free positive 
integer less than or equal to 100, we searched through all the convergents 
with q < 10 100 . _ _ _ 

In this way, we obtain measures for v5, and -^41 — values of n 
within the range considered by Chudnovsky, but not treated by him — as 
well as an improved irrationality measure for \/7. Bennett also found the 
same a and b for these n (along with n = 41 and 57, which we also consider 
here). However, his version of our Theorem 12.11 was not sufficiently strong 
to allow him to obtain effective irrationality measures for n = 41 and 57 
which improve on Liouville's theorem, so these remain as new results here. 

Given the scale of the search, the table is almost certainly complete for 
n < 100. 

The values of a and b listed in Table 1 produced the minimal values of 
k < 2 satisfying the conditions of Theorem 12. II for the given value of n. 
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A key element in translating the sharp result contained in Proposition 
into tight numerical results is a strong bound for 



9(x;k,l) 



E 



log p. 



p — L mod fc;p,primc 

p<x 



Ramare and Rumely [13] provide good bounds. However, due to re- 
cent computational work of Rubinstein [14], we are able to improve these 
bounds considerably for some k. So we present here the following results 
on 6(x; k, I), and the closely-related i/j(x; k, I), for k = 1, 3,4 and 6. 

Theorem 2.3. (a) For 1 < x < 10 12 , 

max max(\e(y)-y\,\ip(y) - y\) < 2.052818^, 

l<y<x 

max max(|%;3,±l) - y/2\ , \if>(y; 3, ±1) - y\) < 1.798158>/£, 

l<y<x 

max max(|%;4,±l) - y/2\ , \if>(y- 4, ±1) - y\) < 1.780719^ and 

l<y<x 

max max(|%;6,±l) -y/2\ , \ip(y;Q,±l) - y\) < 1.798158^. 

l<y<x 

(b) For each (k,l),xo and e given in Table 2, 



6(x;k,l) 



x 



<p(k) 



ip(x; k, I) 



<p(k) 



< ex, 



for x > xq. 





10 5 


10 e 


10 v 


10 8 


10 y 


10 iU 


(1,0) 


0.00474 


0.00168 


0.000525 


0.0001491 


0.0000459 


0.0000186 


(3,1) 


0.00405 


0.00148 


0.000401 


0.0001260 


0.0000371 


0.0000351 


(3,2) 


0.00217 


0.00068 


0.000180 


0.0000428 


0.0000351 


0.0000351 


(4,1) 


0.00494 


0.00169 


0.000471 


0.0001268 


0.0000511 


0.0000511 


(4,3) 


0.00150 


0.00036 


0.000197 


0.0000511 


0.0000511 


0.0000511 


(6,1) 


0.00405 


0.00148 


0.000401 


0.0001260 


0.0000371 


0.0000351 


(6,5) 


0.00217 


0.00068 


0.000180 


0.0000428 


0.0000351 


0.0000351 



Table 2. Analytic epsilons for x > xq 



Only the results for 6(x;3,2) and ip(x;3,2) will be used here, but we 
record the additional inequalities in this theorem for use in ongoing work 
and by other researchers as they improve the current bounds of Ramare 
and Rumely [13] by a factor of approximately 30. 

Unless otherwise noted, all the calculations mentioned in this paper were 
done using the Java programming language (release 1.4.2) running on an 
IBM-compatible computer with an Intel P4 CPU running at 1.8 GHz with 
256 MB of memory. Source code for all programs can be provided upon 
request. Many of these computations were also checked by hand, using 
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MAPLE, PARI/GP and UBASIC. No discrepancies beyond round-off error 
were found. 

3. Arithmetic Properties of Hypergeometric Polynomials 

We use 2F1 (a, 6; c; z) to denote the hypergeometric function 



2-Fi(a,6; c;z) 



1 + 



Aa(a + 1) 



(a + k - 1)6(6 + 1) • • • (6 + k - 1) 



k=l 



c(c+l)-"(c + k— l)k\ 



For our purposes here, we are interested in the following functions, which 
we define for all positive integers m, n and r with (m,n) = 1. Let 

z r iF\ (— r, — r — m/n; 1 — m/n; z^ 1 ) , 
iF\ (— r, — r — m/n; 1 — m/n; z) and 
(m/n) • • • (r + m/n) 



X-m,n,r (-2*) 
^m,n,r (-2) 



i? 



m,n,r 



» 



-2-F1 (r + 1 - m/n, r + 1; 2r + 2; 1 - z) . 



(r + 1) • • • (2r + 1) 

This differs from (4.3) of [5j where the expressions for X r (z) and l^(z) 
have been switched. The same change must be made in (4.4) of [5j too. 

Notations. We let D mnr denote the smallest positive integer such that 
D m ,n,rYm,n,r(z) has rational integer coefficients. 

To simplify the notation in the case of m = 1 and n = 3, which is of 
particular interest in this paper, we let X r (z),Y r (z), R r (z) and D r denote 
Xi ) 3,r(z),Y 1:3)r (z),R lt3)r (z) and -Di, 3 ,r, respectively. 

We will use v p (r) to denote the largest power of a prime p which divides 
into the rational number r. 

Finally, we let [ - J denote the floor function which maps a real number to 
the greatest integer less than that number. 

We first need a refined version of Chudnovsky's Lemma 4.5 in order to 
establish our criterion for the prime divisors of -Dm,n,r- 

Lemma 3.1. Suppose that m,n,p,u and v are integers with < m < n 
and (m,n) = (p,n) = 1. For each positive integer, i, define the integer 
1 < ki < p 1 by kin = m mod p 1 . Then 



m 



E 

i=l 

oo 

E 

i=i 



v - ki 




u — 1 — ki 


pi 




pi 



-u + ki 
pi 



p l 



Remark. It would be more typical to state the above lemma with the 
condition < ki < p 1 rather than 1 < ki < p l . The proof below holds with 
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either condition. However, the above formulation suits our needs in the 
proof of Proposition 13.21 below better. 

Proof. For each positive integer i, we will count the number of j's in u < j < 
v with n j — m = mod p % . That is, with nj — k^n = mod p % . And, since 
(n,p) = 1, with j = k{ modp\ The remainder of the proof is identical to 
Chudnovsky's proof of his Lemma 4.5 [5], upon replacing his p with p l . □ 

Proposition 3.2. Let m,n and r be positive integers with < m < n and 
(m, n) = 1. 

The largest power to which a prime p can divide D m ,n,r is at most the 
number of positive integers i for which there exist a positive integer li sat- 
isfying (k,n) = 1, Up 1 = —m mod n such that 

Up 1 + m ■ (n - li)p % -m- n 

< r modp < . 

n n 

Furthermore, all such i satisfy p % < nr. 

Remark. From the calculations done in the course of this, and other, work 
(see, for example, the notes following Lemmas 13. 3[ I3.4l and l5.lj) . it appears 
that the conditions given in this Proposition provide the exact power to 
which a prime divides D m) n,r- However, I have not been able to prove this. 

Proof. Let a r fo denote the coefficient of z h in Y m ^ r (z) and let p be a prime 
number. From our definition of Y m n r (z) above, we can write 



a r ,h 

where 



r\ a 



r.h 



hi B 



r,h 



B r .h = JJ('" — "' ) and C r ^h = J [ {in + m). 

i=l i=r—h+l 

We first show that if p divides D m , n ,r then (p, n) = 1. 

If p does divide D m ,n,r then p must divide B r ^ for some < h < r. So 
it must divide some number of the form in — m where 1 < % < r. But, if 
p divides such a number and also divides n, then it must also divide m. 
However, our hypothesis that (m, n) = 1 does not allow this and so, if p 
divides D m ^ n ^ then (p, n) = 1. 

Therefore, for any positive integer i, we can find an integer ki with 
1 < ki < p l , (ki,p l ) = 1 and kin = m mod p % . 

As 1 < ki and m < n, we know that < kin — m, and so there must be a 
positive integer li with (li,n) = 1 and kin — m = lip 1 . Furthermore, li < n. 

Returning to our expression for a r ^, we have 

v p (a r>h ) = v p ( ( J J + v p (C r>h ) - v p (B r>h ) . 
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It is well-known that 



oo 

E 

i=l 



r 




h 




r — h 




pi 




pi 




pi 


) 



oo 

E 

1=1 



h — ki 






pi 




pi 



From the first expression in Lemma 13. II with u 

v p (B r> h) 

From the second expression in Lemma[3J]with u 

oo 

Vp {C r)h ) = ^ 
i=l 

Thus, we want to determine when 



1 and v = h, 



-r and v 



-r+h-1, 



( 


r + ki 




r + ki — h 




pi 




pi 



(2) 



r 




h 




r — h 




h — ki 


+ 


ki 


+ 


r + ki 




pi 




pi 




pi 




pi 


pi 


pi 





r + ki — h 
pi 



is negative. 

This will suffice for the purpose of proving this proposition since, as we 
shall show shortly, the expression in ([2]) can never be less than —1. 

We now show that if p % > nr, then the expression in j5J cannot be 
negative. This will establish the last statement in the Proposition. 

Since < h < r < p l for such i, the first three terms in ([2]) are 0. 
Furthermore, the same inequalities for h and r along with the fact that 
ki > show that the sum of the last two terms cannot be negative. 

We saw above that kin — m = lip 1 for a positive integer So it follows 
that kin > p 1 + m > nr > nh. In particular, fc, > h. Furthermore, 
1 < ki < p 1 . Therefore, [(h — ki)/p l \ and [—ki/p l \, are both equal to — 1, 
so the sum of the remaining terms in ([2]) is also zero. 

This establishes the last statement in the Proposition. 

Moreover, if p % > nr + m, then 



r + h < 



P 



n 



m (n — lW + m 
— + — < p\ 

n 



And so the expression in ([2]) is always for such i. We will use this fact 
in the proof of Lemmas 13.31 and 13.41 below. 

For any positive integer i, we can write h and r uniquely as 



h = hup 1 + h i0 and r = rup 1 + r i0 , 



where < hio, rjo < p 1 ■ 
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With this notation, we see that 
r — h 



r 




h 




pi 




pi 





p 



P 



(3) 





h- h 




ki 






P L 






pi 


r + ki 




r 


+ k{- h 


pi 








P 


i 



h a + 1 + 



h i0 - h 



P 



hn + 



no + h 
pi 



and 

rio + h 



Ho 



P 



The first relation holds since 1 < hio, r^o < p % and so [hio/p' l \ = L r io/p*J = 
0. The second relation holding since 1 < ki < p l and so [—ki/p l \ = — 1. 

The last two quantities can only have the values hn or hn + 1, so if 
the expression in ([2]) is to be negative then the first quantity here must be 
zero, since it is never negative, the second must be hn + 1 and the third 
must be ha. This information also substantiates our claim above that the 
expression in ([2]) is always at least —1. 

Since < hio , r ^ < p 1 , the first quantity in (J3j) is zero if and only if 

(4) r i0 > h i0 . 

The second quantity in ([3]) is hn + 1 if and only if 

(5) h i0 > k { . 

Finally, if the last quantity in ([3]) is hn, then [(^o + ki)/p l \ = [(^o + 
ki-h i0 )/p l \. From ©, we find that r i0 + ki-hio < r i0 <p l , sor i0 + ki < p % 
also. Hence 

(6) 



r,n + kj 

< " < 1, 

pi 



the left-hand inequality being strict since ki > 0. 

From ([4]) , we have ki < r^o + h — , while from ([5]) and 



it follows 



that r i0 + ki - h i0 < p l 
find that 



hio < p 1 — ki. Combining these inequalities, we 



k% fio + ki 
pi — pi 



hio < 1 _ h 

pi 



In addition, from ([!]) and ([5]), we know that ki < r^o and from ([6]), 

no < P l ~ h - 1, so 

h + 1 



ki r 

rn + -r < -r < r iX + 1 
pi pi 



P 



Substituting (lip 1 + m) /n = ki into this expression completes the proof 
of the Proposition. □ 

It will be helpful for applications to present a slightly weaker but more 
immediately applicable result on the prime divisors of -D m ,n,r- With that 
in mind, we state the following. 
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ll 



Lemma 3.3. (a) Let r be a positive integer. If p\D m ^ n ^ r , 

\og(nr) 



then 



Vp {Dm,n,r) ^ 



logp 



(b) Ifp is a prime number greater than (nr) 1 / 2 which is a divisor of D m ,n,r, 
then p 2 \D m>n ^ and for some 1 < I < n/2 with (l,n) = 1, Ip = —m mod n, 
and 

m 



nr + m + n nr 
<P< 



^ nA + n-l ' ~ nA + T 

for some non-negative integer A. Moreover, every such prime greater than 
(nr + to) 1 / 2 is a divisor of D m ,n,r- 

Remark. The result in (a) is best possible. E.g., -D2,3,i7 is divisible by 49 
and Ll°§(3 ■ 17)/(log7)J = 2. This example also shows that neither of the 
statements in Lemma [3. 31 holds here (i.e., smaller primes may divide Drn,n,r 
to higher powers than one and p need not lie in the intervals specified by 

©)■ 

Furthermore, 5 divides -D2,3,i0) so the congruence conditions in Lemma [3,3l 
do not hold in general either (since 1 • 5 ^ —2 mod 3). 

Proof, (a) This follows immediately from the last statement in Proposi- 
tion E2J 

(b) Again from the last statement in Proposition l3.2l and our lower bound 
for p, we need only consider i = 1. 

From the inequality on r mod p in Proposition 13.2^ we can write 

Ip + m ^ ^ A ^ i (n — l)p 



(8) 



Ap + 



< r < Ap + 



m 



n 



n n 
for some non-negative A. This provides our upper and lower bounds for p 
in part (b), which suffices to prove the first statement in part (b). 

To prove the second statement, we will show that these primes divide 
the denominator of the leading coefficient of Y r (z). So we let the quantity 
denoted by h in the proof of Proposition 13.21 be r. Using the arguments to 
derive ([3]) in the proof of Proposition 13.21 (|2|) simplifies to 



-1 



m - h 


+ 


r i0 + h 




k% 


= -1 - 


r i0 - ki 




no + h 








+ 


pi 


p % 




pi 




pi 




pi 



where r = r i0 mod p l . 

Therefore, as we saw in the proof of Proposition 13.21 



v p (a rtr ) = 



1=1 





no - ^ 


+ 


no + h 






pi 




pi 



Notice that for i > 2, p 1 > nr + m, so, as we saw in the proof of 
Proposition 13.21 the summands for such i are zero and can be ignored. 



12 



Paul M VOUTIER 



For 4 = 1, from ([8]), (Ip + m)/n < r^o < ((n — l)p — m — n)/n. From 
the relationship between k\ and l\ given in the proof of Proposition 13.21 we 
also have k% = (Ip + m)/n. Therefore, < |_ r io ~~ k\\ , [fio + k\\ < p — 1 
and the summand for i = 1 is -1. Hence v p (a r)T ) = —1, so p divides the 
denominator of a rjr precisely once, completing the proof of the Lemma. □ 

As n gets larger, the structure of the denominator becomes more com- 
plicated and the above is the best that we can do. However, in the case of 
m = 1 and n = 3, 4 or 6, we can obtain a sharper result which will be used 
in this paper. 

Lemma 3.4. Let m = 1 and n = 3, 4 or 6. 

(a) Let r be a positive integer. If ' p\D mjnr then p = n — 1 mod n and 



log (nr) 1 



21ogp 2 

(b) If p is a prime number greater than (nr) 1 / 3 which is a divisor of D m ,n,r 
then p = n — 1 mod n, p 2 )(D m>nr and 

nr + n + 1 nr — 1 

(9) -r——r<P< 



nA + n — 1 ~ ~ nA + 1 
for some non-negative integer A. Moreover, every such prime greater than 
(nr + l) 1 / 2 is a divisor of ' D m ,n,r- 

Remark. The result in (a) is best possible. D42 is divisible by 25 and 
[log(3 • 42)/(21og5) +1/2J = 2. Similarly, D xm2 is divisible by 125 and 
[log(3 • 1042)/(21og 5) + 1/2J = 3. However, it is not true that v p (D r ) > 2 
for all p < (Sr) 1 / 3 (e.g., v 5 {D 43 ) = 1). 

Remark. The second statement in (b) holds for all p > (nr) 1//3 , and this 
is best possible as the example in the previous remark shows, however the 
proof is technical and lengthy. Furthermore, the result here suffices for our 
needs below. 

Proof, (a) We apply Proposition 13.21 As we saw there, (p, n) = 1. For 
these values of n, the only integers less than n and relatively prime to n 
are 1 and n — 1. 

If p = 1 mod n or if p = n — 1 mod n and i is even, then we require 
li = n — 1 mod n to satisfy Zjj> 1 = — 1 mod n. However, with this value of 

Zj , 

(n — 1W + 1 Up 1 + m (n — li)p l — m — n p l — n — 1 
= < r mod p < = 

n n n n 

can never be satisfied. 

If p = n — 1 mod n and i is odd, then we can take l{ = 1. From the last 
statement in Proposition 13.21 p l < nr, so the largest possible i is at most 
log (nr )/ log (p), a fact which completes the proof of part (a). 
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(b) The same argument as for p = n — 1 mod n in the proof of part (a) 
shows that we need only consider i = 1 for p > (rar) 1 / 3 . 

The remainder of the proof is identical to the proof of Lemma f3.3f b). □ 

Lemma 3.5. Let m,n and r be positive integers with (m,n) = 1. Define 

= U^nP 1 ^ and S n , r = H p \ n p V ^ ■ 

(a) Let d be a positive divisor of n. The numerators of the coefficients of 
the polynomials X mritr (l — dz) and Y mnr (l — dz) are divisible by d r . 

(b) The numerators of the coefficients of the polynomials X mnjT (1 — n\i n z) 

— nfi n z) are divisible byn r s n ^ r . 

Proof, (a) This is a variation on part (b) which will prove useful both here 
and elsewhere. Its proof is virtually identical to the proof of part (b). 
(b) This is Proposition 5.1 of [5]. □ 



4. Proof of Theorem [23] 

(a) The bounds for x < 10 12 are determined through direct calculation. 
We coded the Sieve of Eratosthenes in Java and ran it, in segments of 
size 10 8 , to determine all primes less than 10 12 as well as upper and lower 
bounds for 9(x;k,l) and ip(x;k,l) for x < 10 12 . The entire computation 
took approximately 182,000 seconds. 

As Ramare and Rumely note, considerable roundoff error can arise in the 
sum of so many floating point numbers. We handled this issue in a similar 
way to them. We multiply each log by 10 6 , round the resulting number 
down to the greatest integer less than the number as a lower bound and 
round it up to the least integer greater than the number as an upper bound. 
We then sum these integers and store the sums in variables of type long, 
which have a maximum positive value of 2 63 — 1 = 9.233... • 10 18 - a number 
greater than our sums. This is more crude than Ramare and Rumely's 
method, but sufficiently accurate for our needs here. 

In addition to just establishing the desired inequalities, we also compute, 
and have stored, 

(i) our upper and lower bounds for 6*(10 8 i; k, I) and ^(10 8 z; k, I), 

(ii) vr(10 8 i,M), 

. 9(x;k,l) — x/ip(k) . ifj(x;k,l) — x/ip(k) 

(m) mm -= , mm 



ze(10 8 (i-l),10 8 i] \JX ie(10 8 (i-l),10 8 i] 

and 

. . 6(x;k,l) — x/if(k) ^(x; k, I) — x/(p(k) 
fiv) max — , max — , 

xG(10 8 (i-l),10 8 i] V x xe(10 8 (i-l),10 8 i] \JX 

for i = 1,..., 10,000. 

(b) The bounds in part (b) are obtained by applying Theorem 5.1.1 of 
with the L-function zero information calculated by Michael Rubinstein 
We include details of the values used in Table 3, where we round all 
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quantities up by one in the seventh significant decimal (sixth decimal for 
A x , B x , C x , D x for the sake of space). 

Note that for these values of k, there is only one character, X: f° r each 

d. 

For the computation of A x , we followed the advice of Ramare and Rumely 
\13\ p. 414] regarding the evaluation of their K\ and Ki- Using Simp- 
son's rule with an interval size of 0.001 (along with their Lemma 4.2.4), 
we bound from above the integral for K n (z,w) in their equation (4.2.4) for 
u = w . . . 1000. We then apply their Lemma 4.2.3 with w = 1000, which is 
sufficiently large to provide a good upper bound. 

This provides us with an upper bound for e(ip, x, k). 

Using the authors' upper bound for e(6, x, k) on page 420 of [13], we see 
that our results holds for x > xq. 

Proceeding as above, we found agreement with the data that Ramare 
and Rumely present in their Table 1 for k = 1,3 and 4. 



k 


1 


3 


4 


m 


14 


14 


14 


6 


6.289071 • 10 _ ' 


1.256642 • 10" b 


1.798450 • 10~ b 


A(m,5) 


1.082027 • 10 yi 


6.691384 • 10 8b 


4.425147 • 10 84 


R 


2.721552 • 10~ n 


9.085095 • 10- 11 


1.207835 • lO" 10 


e(ip,x,k) 


3.613190 • 10~ b 


7.097148 • 10~ b 


1.001340 • 10" 4 


d 


1 


1 


3 


1 


4 




8000000.365 


4000000.042 


4000000.413 


2800000.0623 


2800000.340 




5.81243 • 10~ 98 


8.94572 • 10" 94 


9.83501 • 10" 94 


1.27527 • 10~ 91 


1.43730 • 10~ 91 




9.09392 • 10" 103 


2.85005 • 10~ 98 


3.04716 • 10" 98 


5.86164 • 10~ 96 


6.37182 • 10~ 96 




7.30396 • 10" 98 


1.13798 • 10" 93 


1.23103 • 10~ 93 


1.63333 • 10~ 91 


1.80646 • 10~ 91 


D x 


1.14495 • 10" 102 


3.63318 • 10" 98 


3.82113 • 10~ 98 


7.52419 • 10" 96 


8.02375 • 10" 96 


E x 


31.414915 


28.3898896 


33.1560902 


26.8928884 


32.8584828 



Table 3. Data for the Proof of Theorem E3 



Since 2.052818 < O.OOOOI862; 1 / 2 for x > 12.2 • 10 9 , the stated inequal- 
ities for 9{x) and ip(x) holds for such x. Using the above sieve code, it 
is straightforward to calculate 9{x) and tp(x) for x < 12.2 • 10 9 . These 
calculations complete the proof of (b) for 9{x) and ip( x )- 

Similarly, 1.798158 < 0.0000351X 1 / 2 for x > 2.7 • 10 9 and a computation 
completes the proof of (b) for k = 3 and 6. 

Finally, 1.780719 < 0.0000511a; 1 / 2 for x > 2.7 • 10 9 and a computation 
completes the proof of (b) for k = 4.. 
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5. Analytic Properties of Hypergeometric Polynomials 

Lemma 5.1. Let r be a positive integer and define N r to be the greatest 
common divisor of the numerators of the coefficients of X r (l — (a — b)x/a), 
where a, b and d are as defined in Theorem \2.1\ 
(a) We have 

1 0.29L>2 



< 



200 r l / Q A r ' 
(b) We have 

< 1.161 • io 39 e °- 91 ^ and (V3) • • • (r + 1/3) 3 dr D r i 40 e o.9iir_ 
N r r\ N r 

Remark. These results are very close to best possible. Chudnovsky [5] 
has shown that D r ~ e ?rv ^ r / 6 = e a9068 -" r as r —* do. 

Remark. We were able to calculate the D r exactly for all r < 2000 (with 
two different methods using both Java and UBASIC 8.8). These actual 
values were equal to the values calculated using Proposition 13.21 This 
strengthens our belief that Proposition 13.21 captures the precise behaviour 
of the prime divisors of D m ,n,r ( a t least for m = 1, n = 3). 

Proof. We will establish both parts of this lemma via computation for r 
up to the point where Theorem 12.31 can be used to prove the lemma for all 
larger r. 

(a) We computed the quantity on the right-hand side for all r < 2000, 
as part of the computation for part (b). We found that its minimum is 
0.00501 . . ., which occurs at r = 13. 

From the second statement in Lemma f3,4( b). we know that if p is a prime 
congruent to 2 mod 3 with (3r + 4)/2 < p < 3r — 1, then p\D r . Since we 
may now assume that r > 2000, we know that (3r + 4)/2 > 3000. 

From Theorem 12.31 and a bit of computation, for x > 3000, we find that 
\6(x; 3, 2) — x/2\ < O.Ollx, so the product of the primes congruent to 2 mod 
3 in that interval is at least e - 71 " 1 - 511 . Therefore, D r /4 r > e 0-014r-i.5H 
Since r 1 / 6 = e^ r ^ 6 < e a0007r for r > 2000, the desired result easily 
follows. 

(b) Here the computation needs to include much larger values of r, so 
we need to proceed more carefully. 

We break the computation into several parts. 

(1) The computation of the factorial and factorial- like product on the 
left-hand side of the second inequality. We shall see below that the product 
of these terms grows quite slowly and they have a simple form, so this 
computation is both easy and fast. 

(2) The computation of 3 dr /N r . From Lemma 13.51 we find that if d = 
then (3, JV r ) = 1, if d = 1 then 3 r \N r and if d = 3/2 then 3 r+V3< - rl ^\N r . For 
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d = 3/2, one can often do better, by directly calculating the numerators of 
the coefficients of X r (l — 3y/3x), by means of equations (5.2)-(5.4) in the 
proof of Chudnovsky's Proposition 5.1 [5]. 

Directly calculating N r is substantially more time-consuming than cal- 
culating 3 r + 1 '3(r!) ) go we always calculate 3 r + 1 '3( r! ) j continue with calcu- 
lating D r and only perform the direct calculation of N r if the size of 
3 3r / 2 D r /3 r+V3( - r ^ warrants it. 

(3) The computation of the contribution to D r from the small primes, 
that is those less than \/3~r , using Proposition 13.21 

To speed up this part of the calculation, and the following parts, the 
primes and their logarithms do not have to be recalculated for each r. 
Instead, we calculate and store the first million primes congruent to 2 mod 
3 (the last one being 32,441,957) and their logarithms before we start the 
calculations for any of the r's. 

(4) The computation of the contribution to D r of all primes from \/3r 
to (3r — l)/(3A(r) + 1) for some non-negative integer A(r), which depends 
only on r. Again, we use Proposition 13.21 as well as the cached primes and 
their logarithms here. 

(5) The computation of the contribution to D r from the remaining larger 
primes. 

From Lemma [3.4( b). we can see that for any non- negative integer A, the 
contribution to D r from the primes satisfying ([9]) changes, as we increment 
r, by at most the addition of the log of one prime, if there is a prime 
congruent to 2 mod 3 between 3(r — 1)/(3A + 1) and 3r/(3^4 + 1), and the 
subtraction of another, if there is a prime congruent to 2 mod 3 between 
3(r — 1) / (3A+2) and 3r/(3^4 + 2). This fact makes it very quick to compute 
the contribution from these intervals for r from the contribution from these 
intervals for r — 1 — much quicker than recomputing them directly. So we 
incorporate this strategy here: for each i < A(r), we store the smallest and 
largest primes in these intervals along with the sum of the logarithms of 
the primes, p = 2 mod 3, in these intervals. 

Again, we use the cached primes and their logarithms for the intervals 
that lie within the cache. 

In this manner, we proceeded to estimate the size of the required quanti- 
ties for all r < 200, 000, 000. This computation took approximately 89,000 
seconds. 

The maximum of 3 dr D r / (N r e°' 911r ) occurs at r = 19, 946 and is less than 
1.161 • 10 39 , while the maximum of (1/3) • • • (r + 1/3)3*" A-/ (]V r e - 911r r!) 
also occurs at r = 19, 946 and is less than 1.176 • 10 40 . 

For r > 200 • 10 6 , we can use the analytic estimates in Theorem 12.31 



Rational approximations ... 17 

From Lemma[331 we know that S^N' 1 < 3 r / 2 -**( rl ). I n addition, r/2- 
ua(r!) < (log r)/(log 3) +0.5 and 

(1/3) •••(r + 1/3) 4 f [ r <te\ 4rV 3 

— 7T eX P l / o ) — n ' 



r! 9 yjx 3x / 

for r > 1, so 

/ s 3 dr , 3 dr (1/3) • • • (r + 1/3) a/ , 

(10) — < 1.8r and — ±-LJ , < 0.8r 4/3 . 

We divide the prime divisors of D r into two sets, according to their size. 
We let D r ^ s denote the contribution to D r from primes less than (3r) x / 3 
and let D r j denote the contribution from the remaining, larger, primes. 

From Lemma 13.4( a). we know that 



Llog(3r)/(21og(p))+l/2j 



p<(3r)!/ 3 

p=2 mod 3 



Now [x/2 + 1/2J < 3/2[x/3j + 1, so 



f 3i/j (\/3r; 3, 2) / 3/ — \] 
D r ,s < exp I * K - ' + 9 (^37; 3, 2 J \ . 

From Theorem l2.3[ and some calculation, we find that 9{x\ 3, 2), tjj(x; 3, 2) < 
0.51x, so 

(11) D r , s < exp (l.28^3r^j . 

From (1101) and (11 ID . we know that 



no\ ^ o 0.000006r OT i Dr,s3 dr (1/3) • • • (r + 1/3) n.000006r 

(12) ^V r _<e and H <e 

for r > 200 • 10 6 . 

We next consider D r ;. 

From Lemma !3.4l fL>) , we see that for any positive integer N satisfying 
3r/(3iV + 2) > (37-) 1 / 3 , we have 

C N N-l *\ 

D r ,i < exp i 0(?r/(3A + 1); 3, 2) - ^ 9(3r/(3A + 2); 3, 2) I . 

[a=o a=o J 

Let i+(x) denote the maximum of 0.5000351 and 9(y; 3, 2)/y for all y > x 
and let t-{x) denote the minimum of 0.4999649 and #(y;3,2)/y for all 
y > x. With the choice N = 200, we can write 

n / fo /^t+(600-10 6 /(3A + l)) ^t_(600-10 6 /(3A + 2))\ 

^- exp rl£ — ^ — £ — ^ — ) 
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since r > 200 • 10 6 . 

With Theorem 12.3( b). we calculate the necessary values of t + {x) and 
i_ (x) and find that 

n ^ O 0.910993r 

for r > 200 • 10 6 . 

Combining this inequality with (|12p yields 

^ <e o.»n r and 3^(l/3)-( r + l/3) 
N r N r r! 

for r > 200 • 10 6 . 

This completes the proof of the lemma. □ 

We now need to define our sequence of approximations to (a/6) 1 / 3 and 
find an upper bound on their size. 

We start with bounds on the size of the polynomials. 

Lemma 5.2. Let m,n and r be positive integers with m < n/2 and let z 
be any real number satisfying < z < 1. Then 

(13) (l + z) r <y m ,„, r (z) < (l + z 1 ' 2 " 

Remark. The upper bound is best possible as can be seen by considering 
z near 0. 

For hypergeometric applications, we are particularly interested in z near 
1, where it appears that the upper bound could be sharpened to 

(2r)\ r(l-m/n) ^ , _ 1/2 ^ 



4T'^-— '-^-(1+z 

r\ T(r + 1 — m/n) 

although we have been unable to prove this. This is an equality for z = 1. 
In the case of m = 1 and n = 3, this extra factor is about 0.8r -1 / 6 . 

Proof. We start by proving the upper bound. 
We can write 

2r 



fc=0 V J 

(r — k + 

m/n) ■ ■ ■ (k — m/n) 



Y r {z) = E a ^ = E ([) ^ ~g + \ + ™' n) n' ' {r ^ m J n) zk 

l — n l — n V / V 



fc=0 fc=0 

We shall show that 



r\ (r- k + l + m/n)---(r + m/n) zk < / 2r \ fc _ 1/2 /2r\ k 



kj (1 — m/n) ■ ■ ■ (k — m/n) \2k—lJ \2k 

This will prove that Y r (z) < (l + z 1 / 2 ) 2r . 
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Since < z < 1, it suffices to show that 
(14) 

fr\ (r — k + 1 + m/n) ■ ■ ■ (r + m/n) ( 2r \ /2r\ 
° fc = \k) (1 - m/n) ■■■(k- m/n) ~ \2k - l) + \2k) = k ' 

We demonstrate this by induction. 

For k = 0, (I14p holds since cto and 6o are both equal to 1. So we can 
assume that (I14D holds for some fc. 
Notice that 

(r — k)(r — k + m/n) 
a ^ = ( k + l-m/n)(k+l) ak aM 
(r - k)(2r - 2k + 1) 

k+1 ~ (ib + ipfe + i) fc ' 



Thus 



afc + i {r — k + m/n) (A; + 1/2) 



(r - fc + 1/2) (fc + 1 - m/n) 6 fc ' 
Since m < n/2, it is apparent that {r — k + m/n)/(r — k + 1/2) < 1 and 
that (k + 1/2) /{k + 1 — m/n) < 1. Since we have assumed that a^/b^ < 1, 
it is also true that 0^+1/6^+1 < 1, which completes the proof of ([14"]) and 
hence the upper bound for Y m ,n,r{z)- 

To establish the lower bound, we again compare coefficients. It is clear 
that ao = (q) and that > (^) for 1 < k < r. Since < 2 < 1, the lower 
bound holds. □ 

Lemma 5.3. Let r be a positive integer, a and b be positive integers with 
b < a. Put 

p r = ^j^-Xrib/a) and q r = ^^-Y r {b/a). 
Then p r and q r are integers with p r q r+ \ 7^ p r+ \q r and 
(15) ^ a + b y< qr < i.iei . 10 39( e 0.9ll 3 -d( o l/2 + b l/^ - 

Ivy* 



Proof. The first assertion is just a combination of our definitions of D r 
and N r along with an application of Lemma 13.51 while the second one is 
equation (16) in Lemma 4 of [2]. 

We now prove the upper bound for q r . 

From Lemma 15.21 

a r Y r (b/a) < (aW + b 1 ' 2 ) . 

The upper bound for q r now follows from Lemma l5.1f b). 
The lower bound for q r is an immediate consequence of the lower bound 
for Y mrijr (z) in Lemma 15.21 □ 
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The next lemma contains the relationship that allows the hypergeometic 
method to provide good sequences of rational approximations. 

Lemma 5.4. For any positive integers m, n and r with (m, n) = 1 and for 
any real number z satisfying < z < 1, 

(16) Z / ^m,fi,v{,^) ^ / m,n,r(^ ; ) — 1) Rm,n,r{.z) • 

Proof. This is (4.2) of [5] with v = m/n. □ 

We next determine how close these approximations are to (a/6) 1 / 3 . 

Lemma 5.5. Let a,b and r be positive integers with b < a. Then 
(17) 



200ag r 



< 



v(o/6) 



1/3 



Pr 



< 



1.176 • 10 40 (a - b) 



e o.9ii 3 -d( a i/2 _ b i/: 



Proof. Using our definitions of p r , q r and R r (z) and the equality expressed 
in Lemma 15.44 we find that 

2r+1 (1/3) • • • (r + 1/3) 



qria/b) 1 / 3 - p r 



a r D r f a\ !/ 3 / a 

N r V6> 



(r + 1) • • • (2r + 1) 



x 2 Fi (r + 2/3, r + 1; 2r + 2; (a - 6)/a) . 

Since {a — b)/a and the coefficients of this hypergeometric function are 
all positive, we have 

2-F1 (r + 2/3, r + 1; 2r + 2; (a — b)/a) > 1. Using the same arguments as in 
the proof of Lemma 15.11 we can also show that 

(1/3) •••(r + 1/3) > 0.29 



(r + l)---(2r + l) ' 4 r rV6' 

Combining these inequalities with the lower bound for q r in Lemma 15.31 
we obtain 

Dr\ 2 0.29(a - 6) 2r (l + 6/a) r a - b 



(18) ^(a/6) 1 / 3 



P7 



> 



4 r r 1 / 6 ag r 

Recall that N r is the greatest common factor of the numerators of the 
coefficients of X r (1 — (a — b)z/a). Since X r {z) is a monic polynomial, N r < 
(a — b) r . The desired lower bound for |g r (a/6) 1//3 — p r \ now follows from 
(|18p and Lemma 15. U fa) , 

To obtain the upper bound, we apply Euler's integral representation for 
the hypergeometric function, we have 



qria/b) 1 / 3 - p r 



D r a T 



i 



2r+l 



(1/3) •••(r + 1/3) /o^ i/a 



7* 



V6 



o 



V(i-trfi-^^ 



-r-2/3 
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Easton (see the proof of his Lemma 8) showed that 



Vp-o'fi-^^"** 



< 



o \ a 

The lemma now follows from a little algebra and Lemma 15.1T b). □ 



6. A Diophantine Lemma 

Finally, we state a lemma which will be used to determine an effective 
irrationality measure from these approximations. 

Lemma 6.1. Let Sgl. Suppose that there exist ko,lo > and E,Q > 1 
such that for allr G N, there are rational integers p r andq r with \q r \ < koQ r 
and \q r 9 — p r \ < l$E~ r satisfying p r q r +i ^ Pr+lQr- Then for any rational 
integers p and q with p/q ^ Pi/qi for any positive integer i and \ q\ > 1/ (2Zo) 
we have 



> , where c = 2ko(2loE) K and k = °^ — . 



Proof. In the proof of Lemma 2.8 of [3], it is clearly noted that this is 
true. The extra Q which appears in the expression for c in the statement 
of Lemma 2.8 of [3] arises only from consideration of the case p/q = Pi/q% 
for some positive integer i. □ 



7. Proof of Theorem EH] 

By the lower bound in Lemma 15.51 we need only prove Theorem 12. II for 
those rational numbers p/q ^ Pi/q% for any positive integer i. 

All that is required is a simple application of Lemma 16.11 using Lem- 
mas [5T3] and [53] to provide the values of ko,lo,E and Q. 

From these last two lemmas, we can choose ko = 1.161 • 10 39 , Zo = 1.176 • 

l 40( a _ 6 )/ 6) E = e -0.911 3 d( a l/2 _ 6 l/2)-2 and Q = e 0.911 .3-d ( fl l/2 + 

Lemma [531 assures us that p r q r +i 7^ Pr+iQr- In addition, Q > e °- 911 3 -1 * 5 
(\/2 + l) 2 > 2.78 > 1 and the condition a > b shows that l Q > 0. If E > 1 
then we can use Lemma 16. 11 

The quantity c in Lemma 16.11 is 



2.322 • 10 39 < 



2.36 -10 40 - 3 d (y/2 + 6 i/2) | K 
e o.9ii5( a i/2 _ 6 i/2) J • 



Under the assumptions that a and b are positive integers with b < a, E > 
1 and k < 2, one can show, by means of calculation and arguments from 
multivariable calculus, that 3 d e" a911 (^/a + y/b)/{b{^fa - y/b)) < 1.822, 
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the maximum occurring for a = 14 and b = 11. So we can simplify the 
expression above, bounding it above by 

2.322- 10 39 (4.3- 10 40 ) K . 

By the lower bound in Lemma l5.5l for the pi/q^s, we know that the c\ in 
Theorem 12.11 will be a constant times a. Furthermore, we know that, a > 5 
is required in order that E > 1 and k < 2. So we can introduce a factor of 
a/5 into our expression for c above, obtaining 

m 40 n 

2.322 -10 39 (4.3 •10 40 ) K < —-(4.3 • lO 40 )" 

V \/2L5 J 

< 10 4 °( K+1 )a, 

since k < 2. 

The condition that E > 1 (so that a/2 < b < a) along with Liouville's 
theorem shows that Theorem 12.11 is also true if k > 2. 

By these estimates and Lemma [6. II we now know that Theorem 12.11 holds 
once \q\ > 1/(21q) > b/ (2.36 • 10 40 (a — 6)). There is a simple argument we 
can use to deal with g's of smaller absolute value. 

If p/q did not satisfy ([T]), then | (a/6) 1 / 3 — p/q | < 1/ (2g 2 ) would certainly 
hold and p/q would be a convergent in the continued fraction expansion of 
(a/b) 1 / 3 . 

Since b < a, it follows that 36 2 / 3 < a 2//3 + (a6) 1 / 3 +6 2 / 3 . As a consequence, 

36 2/3 ( ffl l/3 _ 6 l/3) 

< a — 6, or, more conveniently, 

/a\i/3 _ a 1 / 3 -^/ 3 a- 6 
UJ ~ ~ W~3 < 3b ' 

So we know that the continued fraction expansion of (a/6) 1//3 begins 
[l;x, ...] where x > [36/(a — b)\ . Therefore po = qo = 1 (here Po/qo is 
the 0-th convergent in the continued fraction expansion of (a/6) 1 / 3 ), while 
Qi > L36/(a - 6)J and it is certainly true that q 1 > 6/(2.36 • 10 40 (a - 6)). 

Hence p/q = 1, in which case a/6 > (6 + l)/6 and E > 1 imply that 
(a/6) 1 / 3 - 1 > 1/(46) > l/(8a) and © holds. 

This completes the proof of the Theorem 12.11 



8. Proof of Corollary [231 

To prove Corollary 12.21 we first need to obtain obtain a lower bound 
for \^/n — p/q\ from the irrationality measure we have for the appropriate 
\Ja/b. There are two different ways of doing this. 
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(i) If v/ a/b is of the form s^/n/t then, from Theorem 12.11 we obtain 



sp 

17i 



> 



and, as a consequence, 



n 



> 



Ci\tq\ K - 



1 



scit K |g| K+1 

Let us look at the case of n = 2 to see how we proceed here. We have 
a = 128, b = 125, s = 4 and t = 5. From Theorem [231 we have c x = 2 ■ 10 97 
and k = 1.4321, so 

io-" 



■iT2- v - 
q 



> 



12.4321 



by the above reasoning. 

We wrote a program in Java to calculate the first 500, 000 partial frac- 
tions and to bound from below the denominators of the first 500, 000 conver- 
gents in the continued-fraction expansion of y/2. For this, we used the algo- 
rithm described by Lang and Trotter [TU] which uses only integer-arithmetic 
and does not require any truncated approximations to \[2. 

The denominator of the 500, 000-th convergent is greater than io 257 ' 000 
and it is easy to verify that 

10-" 0.25 



12.4321 



> 



12.4325 



IQI " l<?l' 

for all q whose absolute value is larger than that. Thus, it only remains to 
check that the desired inequality is satisfied for all q whose absolute value 
is at most the denominator of the 500, 000-th convergent. 
Rather than actually checking directly to see if 

0.25 



^2- 



V 



> 



\q 



2.4325 



held for all these the convergents of \/2, we simply looked at the partial 
fractions in the following way. 

From the theory of continued-fractions, one can show that 



1 



< 



a 



P, 



(o*+i + 2) 

where Oj+i is the i + 1-st partial fraction in the continued- fraction expansion 
of a while pi/qt is the i-th convergent. 

As we see in Table 4, the largest partial fraction found for \[2 was 
&484708 = 4,156,269. Therefore, the corollary holds for |g| > 9 • 10 13 > 
((4156269 + 2)/4)( 1 /°- 4325 ). Now a direct check among the smaller conver- 
gents completes the proof of the corollary for n = 2 (the constant C2 = 0.25 
arises here). 
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We proceeded in the same way for n = 9, 10, 18, 19, 20, 22, 25, 28, 30, 
36, 43, 44, 49, 57, 65, 66, 67, 68, 70, 76, 83 and 84. 

(ii) The other possibility is that {/a/6 is of the form s/{t^/n). In this 
case, we use the fact that \\/x — l/y\ = \ (x — y)/{xy)\ and find that 

.3/ 



3/ n 



n 

> 



sc\\p\t K \q\ R 

We can assume that \ ^/n — p/q\ < 0.5 and so 

3/ 



3/ n- 



P 



P 



"n 

> 



scit K {^i+ l/2) K |p| K + 1 ' 

We then proceed in the same way as in the previous case. 

It is in this way that we prove the Corollary for n = 3, 4, 5, 6, 7, 11, 12, 
13, 15, 17, 26, 31, 37, 39, 41, 42, 50, 52, 58, 60, 61, 62, 63, 78, 90, 91, 98 
and 100. 



[i 

[2 
[3; 

[5 

[7; 

[8 

[9 

[10: 

[11 

[12 

[13 

[Id 
[15. 
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n 


a 


b 




i 






2 


2-4 3 


5 3 


484, 


708 


4,156, 


269 


3 


3^ 


2 A 


13, 


628 


738, 


358 


4 


2-4 3 


5 3 


485, 


529 


8,312, 


539 


5 


239645788 3 


5 • 140145707 3 


266, 


405 


3, 494, 


436 


6 


467 3 


6 • 257 a 


238, 


114 


466, 


540 


7 


44 3 


7- 23 d 


274, 


789 


12,013, 


483 


9 


9 


2 y 


97, 


298 


1,063, 


588 


10 


5- 13 a 


4. 143 


371, 


703 


1,097, 


381 


11 


25022 3 


11 • 11251 a 


217, 


358 


1,352, 


125 


12 


9- 29 a 


4-38 3 


34, 


767 


1,185, 


798 


13 


57 a 


13 • 37 3 


55, 


205 


1,406, 


955 


15 


5 2 


3- 2 A 


245, 


733 


1,571, 


507 


17 


18 a 


17- 7 a 


169, 


765 


1,536, 


142 


18 


9 -29 s 


4-38 3 


300, 


238 


3,143, 


844 


19 


19 • 3 a 


8 3 


138, 


226 


521, 


398 


20 


20 • 7 a 


19 a 


72, 


509 


1,840, 


473 


22 


11 • 5 a 


4. 73 


232, 


141 


595, 


645 


25 


239645788 3 


5 • 140145707 a 


20, 


862 


2,449, 


303 


26 


3 3 


26 


252, 


311 


1,722, 


109 


28 


28 


3 3 


275, 


575 


1,654, 


773 


30 


10 


9 


228, 


793 


197, 


558 


31 


22 3 


31 • 7 3 


205, 


544 


1,643, 


436 


36 


467 3 


6 • 257 a 


238, 


549 


2,799, 


247 


37 


10 a 


37 -3 a 


494, 


731 


6,591, 


064 


39 


39 z • 2 6 


23 d 


309, 


275 


483, 


161 


41 


100 3 


41 • 29 3 


321, 


697 


417,960, 


093 


42 


49 


6- 2 A 


408, 


968 


409, 


489 


43 


43 • 2 3 


7 3 


227, 


706 


1,359, 


766 


44 


44 -2 3 


7 3 


260, 


709 


370, 


994 


49 


44 3 


7- 23 3 


273, 


736 


1,716, 


211 


50 


20 • 7 3 


19 a 


54, 


577 


2,055, 


429 


52 


2 • 2 3 


13 


379, 


989 


3,958, 


641 


57 


57 • 33 d 


127 d 


110, 


601 


847, 


651 


58 


4-2 3 


29 


172, 


932 


139, 


963 


60 


2 • 2 3 


15 


44, 


247 


461, 


876 


61 


4 3 


61 


76, 


517 


3,405, 


348 


62 


4-2 3 


31 


400, 


816 


330, 


326 


63 


4 3 


63 


168, 


229 


2,664, 


200 


65 


65 


43 


183, 


363 


16,950, 


688 


66 


33 


4- 2 A 


179, 


933 


589, 


781 


67 


67 


4 3 


419, 


845 


937, 


766 


68 


17 


2 • 2 4 


121, 


095 


1,059, 


335 


70 


35 


4- 2 4 


376, 


116 


582, 


245 


76 


19 • llll 3 


2 ■ 2353 3 


300, 


013 


575, 


574 


78 


47 3 


78 • ll 3 


421, 


553 


1,145, 


724 


83 


83 • 58 3 


253 a 


431, 


244 


434, 


543 


84 


84 • 33856 3 


148273 3 


236, 


330 


5,018, 


560 


90 


3 • 3 3 


10 • 2 :i 


43, 


615 


314, 


175 


91 


9 3 


91 • 2 3 


123, 


567 


416, 


579 


98 


28 


3 3 


274, 


960 


23, 166, 


836 


100 


5 • 13 3 


4. 143 


336, 


362 


1,383, 


591 



Table 4. Data for the Proof of Corollary 12.2 



